GEOMETRY OF WHIPS AND CHAINS 
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<£ ; 1. Introduction 

The purpose of this note is to explore the geometry of the inextensible 

string (or whip) which is fixed at one end and free at the other end. In 

v— / | the absence of gravity, the inextensible string is a geodesic motion in 

H ' the space of curves in IR 2 parametrized by length and passing through 

,£3 ■ a fixed point. This is a submanifold of the space of all curves in M 2 , 

and the Lagrange multiplier from the constraint is the tension in the 
string. Hence the string satisfies a wave equation, where the tension 
is generated by the velocity implicitly by solving a one-dimensional 
^h ■ Laplace equation. 

The situation seems closely analogous to the motion of an incom- 
pressible fluid, where the configuration space is the space of volume- 
lO ■ preserving diffeomorphisms, a submanifold of the space of all diffeomor- 

phisms, and the Lagrange multiplier is the pressure. This is also a geo- 
desic motion, but the Riemannian geometry is still not well-understood. 
qq | It is our hope that understanding the one- dimensional situation better 

will help in understanding the two- and three-dimensional versions. 
>• . Part of the problem with the inextensible string is that the geo- 

desic equation is not an ordinary differential equation on an infinite- 
dimensional manifold (unlike the situation for incompressible fluids). 
The derivative loss cannot be avoided, and thus one must treat it as 
a hyperbolic partial differential equation. However, since the tension 
must be zero at the free end, the equation degenerates at the free end, 
and the standard approach to symmetric hyperbolic equations must 
be modified to deal with this. Although these issues are not discussed 
here, we will study them in a forthcoming paper. 

The other serious issue is that since the tension is not given by a 
positive function of the speed, position, and time, but rather given 
implicitly by the solution of an ordinary differential equation, it is pos- 
sible for tension to become negative. This corresponds to the differen- 
tial equation changing type from hyperbolic to elliptic, and drastically 
complicates the problem. We show that in the absence of gravity, the 
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tension is always nonnegative, while with gravity, the tension is non- 
negative as long as the string hangs below the fixed point. 

To better understand the necessary conditions on weak solutions 
(which has remained controversial in recent papers on falling strings 
with kinks), we consider the corresponding finite-dimensional model: 
a system of iV particles linked by rigid massless rods. In the absence 
of gravity, this is also a geodesic motion on an iV-torus with a certain 
Riemannian metric. The geodesic equation ends up being a sort of dis- 
cretization of the wave equation for the inextensible string, which still 
conserves energy. Furthermore, the tension in the chain converges to 
the tension in the corresponding whip as iV — > oo. As a result, the Rie- 
mannian curvature of T^ converges to the Riemannian curvature of the 
space of curves, and one can therefore understand the geometry of the 
infinite-dimensional space quite well by studying the finite-dimensional 
approximation. 

This suggests that one may be able to use the same approach to 
study incompressible fluids. Namely, take a collection of particles on 
a torus T 2 in a grid, and constrain them so that the areas of quadri- 
laterals are constant. (One could consider constraining triangle areas, 
but this seems too rigid; it does not appear possible to approximate 
a general area-preserving diffeomorphisms by such maps.) The curva- 
ture of the finite-dimensional manifold should approach that of the full 
area-preserving diffeomorphism group, and hence one expects one could 
understand the geometry of fluid mechanics well by a finite-dimensional 
approximation. 

2. Discrete string 

Let us compute the geometry of a discrete string (a chain), which 
we model by n + 1 point masses joined in M 2 by rigid rods of length - 
whose mass is negligible. We assume the th point is fixed at the origin 
while the n th point is free. 

The unconstrained Lagrangian is 

1 n n 

£= 2 5^l 2 ~5Z^ Xi ' e2 ^ 

i=\ i=l 

In addition the constraints are given by 

,2 1 

Xj -Xj_i = — , 1 < % < n. 

Therefore the equations of motion are 
(1) ^ = -ge 2 + Aj+i(x m - Xj) + Ki^k-i ~ *i) 
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for 1 < % < n, where we are fixing xo = and A n+ i = to make the 
equations satisfied also at the endpoints. 

The constraint equations determine the Lagrange multipliers A, which 
are obviously just the discrete tensions. We get 



Vi-ViJ 2 = -A,; 



2 

l\Xj— Xj_i,Xj + i— Xj) + — Aj — Aj_i(Xi— Xj_i,Xj_i — Xi_2/ 

for 1 < % < n (again using A n+ i = 0), while for i = we get 

|vi| 2 - g(x 1 ,e 2 ) = -A 2 (xi,x 2 -xi) + — Ai. 

We can simplify these A-equations somewhat if we change variables. 
For each 1 < i < n, let x.; = Xj_i + ^(cos^sin^), where 9i G S 1 . 
Then 

(xj - Xi_i, x i+ i - Xi) = — cos (6 i+1 - &i) 



while 



Wi - v z _i 



I 2 = — \9-\ 2 
ri 2 



So the constraint equation for A becomes 

(2) - cos (0 i+1 - ^)A< + i + 2A, - cos & - 9 i . 1 )X i . l = \9\\ 2 
for 2 < i < n and 

(3) -cos(# 2 -#i)A 2 + Ai = \e x \ 2 -ngsvaOx. 
For the evolution equation, we have 

(4) Bi = X i+ i sin (6 i+1 - t ) - A<_i sin (^ - 0^) 
for 2 < i < n and 

(5) #i = A 2 sin (#2 — #i) — ng cos 9i. 

Our first question is whether the tension is positive. We can write 
equations ([3]) and (|2]) as 



^ MyAj = |^| 2 - ngdn sin9 1 , 

7 = 1 



where the matrix M is 

i -cos(e a -^i) o 

-cos(e 2 -6»i) 2 -cos (6» 3 -0 2 ) 

-cos(0 3 -0 2 ) 2 
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If we let <2j = cos ( 
elimination yields 



-i — 6i) for 1 < i < n— 1, then standard Gaussian 



/ 6i -ai 

fe — 0,2 
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a 2 


b 1 b 2 


H 



h 1 rr n— ! a * rr n_1 a * rr n_1 a » 1 1 

°n I lli=l b" lli=2 67 lli=3 tj —J-/ 



where the sequence bi is defined by b\ = 1, 6 i+ i 
for every z, we see that every bi satisfies 1 < fej < 2. 
Reducing further, we obtain the inverse matrix 

n _. m— 1 m— 1 

(6> «*= e ^n?rir- 



2 — X-. Since laj < 1 

^z 



m=maxtj fc=« t=j 



Thus the tensions are 



\ i = M il {\0 l \ 2 -ngsin0 l ) + ^M i] \0 3 \ 2 . 

3=2 

If every a^ > 0, then every component of M _1 is nonnegative. So as 
long as \0i+i — 0i\ — n f° r every 1 < % < n — 1 and sin^ < 0, we will 
have all tensions \ > 0, regardless of \0i\. If any a^ < 0, then some 
M u < 0, and thus if 0k = Sjk, then Aj < 0. Additionally, we always 
have M 11 > 0, so that if sin#i > 0, then when all 6\ = 0, we have 
Ai = -ngM 11 sin 1 . 

Let us now look at the curvature; we will assume for the moment 
that g = to get a purely geometric problem. We consider M n = (S* 1 )" - 
as a submanifold of R 2n , with induced Riemannian metric given by 
the kinetic energy formula above. We know that the equation of a 
constrained geodesic is always 



d 2 x k 
dt 2 



dk = Bi. 



1 F). I f) 

dt l) dt J 



d k . 



For any vector u G TgM n , we have (v,k — Uk-i,Xk — Xk-i) = 0, for every 
1 < k < n. So Uk — Uk-i = -^r]k(— sin 0k, cosOk) for some rjk G R. In 
terms of the 77's, we can write (using (J4])) 

(7) B fc (tt, tt) = \{U, M)fc+l(x fc+ i - X fe ) + X(u, M)fc(x fc _i - x fc ), 

for 1 < k < n, where 

n 

(8) A(u,u), = ]£m<%| 2 . 

3=1 
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Since the ambient space is flat, the curvature is, by the Gauss- 
Co dazzi formula, 

(R(u, v)v, u) = (B(u,u),B(v,v)) — (B(u,v),B(u,v)). 

To simplify this, we first compute 

1 n 
(B(u, u), B(v, v)) = — ^2 A ( M ' w)fc+iA(u, v) k+l + X(u, u) k X(v, v) k 

fc=i 

- X(u, u) k+1 X(v, v) k cos (9 k+1 - 9 k ) 

- X(u, u) k X{v, v) k+1 cos (0 k+1 - 9 k ) 

= —X(v, v)i (X(u, u)i - X(u, u) 2 cos (9 2 - 9i)\ 

1 n 
+ —z V^ 2A(m, u)kX(v, v) k 

- X(u, u)k-iX(v, v)k cos (9 k - 9 k -i) 

- X(u, u) k+1 X(v, v) k cos (6 k+1 - 9 k ) 

1 n n 

= -2-J2Yl M *t(0)Hu, u)iX(v, v) k 
fc=l l=l 

t=i j=i 

Similarly we have 

.. n n 

(B(u,v),B(u,v)} = -J2J2 MiJ ( e ^iVjtr 
n »=i i=i 

Therefore the curvature is 

1 n 
(9) (R(u, v)v, u) = -J2 Ml3 ($) fa?£ ~ V&V&] 






1 J2MH9)[U J - V ^ l \ 



2n 2 

Thus the curvature is nonnegative in every section if and only if ev- 
ery M l i(9) > 0, which also happens if and only if every tension is 
nonnegative regardless of the initial velocity. 

To get a bound on the curvature, we want a bound on the components 
M !J . We see from equation ([6]) that M y < max(n + 1 — i, n + 1 — j). 
The extreme case is when all 0's are the same (so that Oj = 1 for 
every i), which corresponds to a straight string. Then we have M l i = 
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max(n + 1 — i,n + 1 — j). Therefore the sectional curvature satisfies 
< K < n. 

3. Smooth string 

Now let us consider the continuum model of an smooth inextensi- 
ble string. The configuration space is the set Cl,o of all C°° maps 
77: [0,1] -> M 2 such that 77(0) = and (rf(s),r/(s)) = 1. This is a 
closed submanifold of C = {r] G C°°([0, l],M d ) 1 77(0) = 0}. The tan- 
gent spaces are 

T V C = {X G C°°([0, 1],K 2 ) I X(s) G T„ (s) M 2 Vs, X(0) = 0} 

and 

T V C L , = {XE T V C I (X'(s),r ] '(s)) = Vs}. 
The Riemannian metric is defined by 

((X,Y))= ! (X(s),Y(s)) v{s) ds. 
Jo 

In this metric, the orthogonal space to T v Cl,o is 

(T,C L)0 ) ± = {Y G C°°([0,1],R 2 ) j F( S ) = l(er( S )^( S )) 

for some a G C°°([0, 1],R) with or(l) = 
The orthogonal projection of a vector X G T^Co should be 

Pl(X) = X - ±(<xrf), 

as 

where a satisfies 

(10) a"(s) - (rj"(sW(s))a(s) = (X' \s) , rj' \s)) , a(l) = 0, 

but we have a problem with the boundary condition at s = 0. We want 
Pe, to map into T^C^o, but 

P L (X)(0) = X(0) - a'(0W(0) - a(0)rj"(0) = -a' (0)rj' (0) - a(0)ri"(0). 

Now r)"(0) is orthogonal to r)'(0), so if rj"(0) 7^ 0, we need both <r(0) = 
and cr'(0) = 0, which will generally be impossible. So we do not have 
a nice orthogonal projection. There are two ways to get around this. 

(1) Consider only curves with r/'(0) = 0. This forces X"(0) = as 
well, and then we have to worry about Pl(X) also satisfying 
(P L (X))"(0) = 0. This forces 77* (0) = and X iv (0) = 0, and 
so on... We end up being forced to require that all derivatives 
of rj exist at s = (which is why I assumed 77 G C°°([0, 1], M 2 )) 
and that all even derivatives of 77 vanish at s = 0. Alternatively, 
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r] extends to a C°° curve on [— 1, 1] such that T)(—s) = —fj{s). 
This observation leads to the second idea. 
(2) We could assume from the start that r\ is the restriction of an 
odd curve on [—1,1]. Then the tangent space will automati- 
cally require X odd. We then have two boundary conditions 
for (dO]): o-(-l) = and <r(l) = 0. Finally, the solution of dTOj) 
will automatically be even in s by symmetry, and then we can 
just restrict to [0,1]. The advantage of this approach is that 
we don't need to assume any more smoothness than is abso- 
lutely necessary for (1101) to make sense. This is of course the 
same kind of trick that's used to get simple solutions to the 
constant-coefficient linear wave equation with one free and one 
fixed boundary. 

Either of these two tricks will force cr'(0) = 0, so we may as well assume 
this as our other boundary condition. 

These difficulties do not arise if both ends of the curve are free. 
If both ends of the curve are fixed, we need oddness through both 
endpoints, again as with the simplest wave equation. 

Now, the first thing we do with the orthogonal projection is to com- 
pute the second fundamental form B: T v Cl,o x T v Cl,o - > (TjjCl.o) 
defined by B(X, Y) = V X Y - P L (V X Y). We obtain 

(11) B(X,Y) = ±{a XY (sW(s)) 

where 

(12) a' XY {s) - K 2 (s)a XY (s) = -(X>(s),Y'(s)). 
This formula implies the geodesic equation 

d 2 rj ( di] dr]^ 



dt 2 \dt' dt 

or more explicitly, 

r)tt(s,t) = dg(a X x(s,t)r] a (s,t)), 

d 2 cr xx (s,t) - \r] ss (s,t)\ 2 o- xx (s,t) = -\r] st (s,t)\ 2 , 

with boundary conditions as discussed. 

If we want to incorporate gravity, then we have 

(13) r]tt{s, t) = d s (cr(s, t)r] s (s, £)) - ge 2 

(14) o- ss (s,t) - \rj ss (s,t)\ 2 cr(s,t) = -\r]st(s,t)\ 2 , 

with boundary condition cr(l, 0) = 0. To determine the boundary 
condition at s — 1, we assume r/ ss (0, t) = as before. If we want 
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T) tt (0, t) = 0, we must have a s (0, t)r) s (0, t) = ge 2 , which cannot be sat- 
isfied unless e 2 || r] s (0,t). Instead, motivated in part by the finite case 
above, we set <r s (0,t) = g(e 2 ,77 s (0,£)), which is equivalent to forcing 
the acceleration parallel to the curve to be zero. 

If, as before, we use the fact that \r) s \ = 1 to write rj s (s,t) = 
(cos9(s,t),sm9(s,t)), then differentiating the first of equations ffl3l) 
with respect to s yields 

(15) 9 tt (s, t) = 2a s (s, t)9 s (s, t) + a(s, t)9 ss (s, t) 

(16) (T ss (s, t) - 9 s (s, t) 2 a(s, t) = -9 t (s, tf 

with boundary conditions <r s (0,t) = g cos 9(0, t) and cr(l,t) = 0. 

Now we want to compute the curvature; again we assume g = to get 
a purely geometric formula. Formula (TTT|) gives, via the Gauss-Codazzi 
equations, the curvature formula: 

({R(X, Y)Y, X)) = ({B(X, X), B(Y, Y))) - ({B(X, Y), B(X, Y))) 

1 / d d 

{<Txx(s)r)'(s)),—(aYY(sW(s)) 



i, \ds ds 



-^((TxY(s)r]'(s)),—(a X Y(s)7]'(s)) ) ds 



= / cxx(s) ( - °yy( s ) + K 2 {s)a YY {s)) ds 
Jo 

- / ffxy(s) ( - 4y( s ) + K 2 (s)(Jxy(s)) ds. 
Jo 

If we denote the Green function of equation (|12l) by G(s, q), so that 

a XY (s)= [ G(s,q){X'(q),Y'(q))dq, 
Jo 



then the curvature formula becomes 

-1 r i 



{(R(X 1 Y)Y,X))=J J G(s,q)(\X'( q )\ 2 \Y'(s) 

-(X'(q),Y'(q)){X'( S ),Y'(s)))dsdq 



1 rl rl 



G(s,q) J2 {X'{s)Y;{q) - X'{q)Yl{s)) 2 dsdq, 



2 - /l! • / " ij=i 



using the symmetry of the Green function. Thus if G(s, q) > for all 
s and q, we get nonnegative curvature. 
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To prove positivity of G(s, q), we first compute: 

G{q,q) = \ G(s,q) S(s — q) ds 
Jo 

-i 

G(s, q) (G ss (s, q) - K 2 (s)G(s, q)) ds 



f G s (s,q) 2 + K 2 (s)G(s,q) 2 ds > 0. 
Jo 



So G(q, q) > 0. We also know that G s (s = q + , q) —G s (s = q~, q) = — 1. 
Since the boundary conditions are G s (0,q) = and G(l,q) = 0, we 
can say the following. 

(1) If G(0, q) < 0, then G(s, q) < for s near 0. As long as 
G(s,q) < 0, we have G ss (s,q) < 0, which implies G s is non- 
increasing. Since G s (0, q) = 0, this means G s (s,q) < as long 
as G(s, q) < 0. Thus there is no way G(s, q) can increase to 
or any positive number on < s < q. This makes G(q, q) > 
impossible, a contradiction. Therefore G(0, q) > 0, and further- 
more G(s, q) > for all s G (0, q). 

(2) If G s (l,q) > then G(s,q) < for s close to 1. As long as 
G(s, q) < 0, we have G ss (s, q) < 0, and so G s (s, q) is nonincreas- 
ing. So G s (s,q) > as long as G(s,q) < 0, so G cannot have 
a turning point on q < s < 1. This again makes G(q,q) > 
impossible, implying a contradiction. So G s (l,q) < 0, which 
implies G(s, q) > for all s G (q, 1). 

This proves that the curvature is nonnegative as long as rj is smooth. 
In fact the sectional curvature must be strictly positive in all nontrivial 
two-planes, since it can only vanish if there is a c such that X'(s) = 
cY'(s) for every s. Since X(0) = Y(0), we must have X(s) = cY(s) for 
all s, which means X and Y are linearly dependent in T v Cl^. 

Polarizing the formula for curvature, we obtain 

((R(Y, X)X, W)) = J X j X G(s, q) (\X\q)\ 2 (Y'(s), W\s)) 

-(X'(q),Y'(q))(X'(s),W'(s)))dsdq 
1 r l 



JO 



•V'«//): 2 {§- s {G(s,q)Y'(s)),W(s: 



+ (X'(q),Y'(q)) (^- s {G(s,q)X'(s)),W(s))dsdq. 
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Therefore the curvature operator is 



Y " p ^.lTs 



d 



G( S , q )(-\X'( q )\ 2 Y'( S ) + (X'( q ),Y'( q ))X'(s))dq . 

Even for smooth X, this operator is not bounded in L 2 . (By contrast, 
the curvature operator in the volumorphism group for a fixed C 1 vector 
field X is bounded in L 2 and in any other Sobolev space.) 

4. Convergence 

We don't have a good existence and uniqueness result for the continu- 
ous equation; the closest we have is Reeken [R2J , who proved short-time 
existence for a chain in M 3 in something like the Sobolev space H 17 with 
initial conditions close to a chain hanging straight down (for technical 
reasons, Reeken assumes the chain is hanging from (0, 0, oo) with non- 
constant gravity). This result is already quite difficult. It is conceivable 
that we could get a simpler result using the two-dimensional equation 
derived above. Of course, the finite-dimensional approximation is the 
flow of a vector field on a compact manifold (a bounded subset of the 
tangent space T^S 1 )" 1 , since the kinetic energy is bounded). Thus the 
finite-dimensional approximation has a unique solution existing for all 
time and for any n. We can thus ask whether solutions of the PDE can 
be constructed as limits of the ODE as n — >• oo. 

It is easy to see that if we knew a smooth solution of the PDE existed, 
then the discrete model would converge to it. Let us write 

6 n : i -,-,..., li xt^ S 1 andff„: i-,~...,llxR->R. 
yn n ) yn n J 

We will set 6 n (%,t) = 9 k (t) and <r n (J,t) = ±\ k (t) for 1 < k < n. 
Then equations (J2J) and (jlj) become 

(17) p n (^t) =n 2 a n (^t) S in[6 n (^,t) -e n (^t)} 

-n 2 a n (^,t)sm[e n ^,t)-9 n ( k -^,t)] 
and 

(18) -cos[6 n ( k -^,t) -6 n (^t)]a n (^t) +2a n (^t) 

-cos [e n ^t)-e n (^t)]a n (^t) = i Q^,t)) 2 

The operators appearing in equations (TT71) and (TTB1 are symmet- 
ric discretizations of derivative operators; if a and 9 are sufficiently 
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smooth, we can write 

n 2 a(x + i,t) sin [9(x + J,t) - 9(x,t)] 

-n 2 a(x-±,t) sin [9(x,t) -9(x-±,t)] 

= a(x)6"(x) + 2a'(x)9'(x) + 0(4,) 

and 

- n 2 cos [9(x + ±,t) -0(x,t)]<r(x+±,t) + 2n 2 a(x,t) 
-n 2 cos [9( x ,t)-6(x-±,t)]<T(x-±,t) 

= 9'(x) 2 a(x)-a"(x) + 0{^) 

where the constant in 0(\) depends on the C 4 norms of both 9 and 
a. 

Therefore, we can say that if a and 9 form a sufficiently solution of 
( TT5j) . then the motion of the discrete chain with n elements converges 
to the motion of the smooth chain as n — > oo. We have convergence not 
only of the position and velocity, but also of the acceleration and the 
tension, which is not necessarily typical (for example, a system with 
strong returning force to a submanifold has position and velocity con- 
verging to the geodesic motion on the submanifold, but the acceleration 
does not converge). 

Because the convergence is so strong, we can speculate about conver- 
gence of the motion in the nonsmooth case. Suppose the continuous 
string has a kink in it. There has been some serious analysis of the 
possible jump conditions one should impose on the equations in this 
case. See for example Reeken |Rlj . who used energy conservation and 
momentum conservation to classify possible motions of a kink. Also 
D. Serre used an approach similar to yours, relaxing the constraint to 
|^'(s)| < 1 and requiring that the tension always be nonnegative, and 
derived jump conditions. But this only sets up the differential equa- 
tions, and there's no guarantee that they have solutions. (Serre found 
blowup phenomena in a couple of examples.) 

Indeed, since we expect the motion of n constrained particles to con- 
verge to the motion of a chain, the fact that the particles can (ideally) 
support negative tension in certain motions (when there is an acute- 
angle kink or when the chain starts with the first link above the fixed 
point) seems to suggest that the motion of the continuous string should 
also support negative tension. Of course this destroys the actual equa- 
tions, which suggests that the motion of n constrained particles has no 
limit in some cases. 
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We will look at the limit as n — > oo by examining the inverse ma- 
trix M l i given by §§§. We expect the components of this matrix to 

converge to the Green function as -M* J « G[ ± ,-). To understand 

what happens to M % \ we first look at the sequence hi appearing in its 
definition. We have 

a 2 
b i+1 = 2 - -y-, where a t = cos (8 i+1 -6i). 
hi 

Now if the discrete curve approaches something smooth, then we will 
have 

&i+l — 9% ~ —K{, 

n 

where Ki = 9( l — ), the approximate curvature. Thus as long as the 
curve is smooth, we will have at ~ 1 and thus also 6j « 1. 

Numerically experimenting, we find that the sequence 6« is approxi- 
mated very well by 

where / satisfies the differential equation 

(20) /'( S ) = k 2 ( S )-/ 2 ( S ), 

with initial condition /(0) = 0. This equation of course comes from 
the Riccati trick f(s) = y'(s)/y(s), where y"{s) — K 2 (s)y(s) = 0. 

Now suppose there is a kink in the curve, say of angle a at position 
s with s ~ k/n for some integer k. Then k+ i —6 k = a, and so all aj's 
are close to 1 except for a k = cos a. We then have b k +i being rather 
far from 1, which means for equation ( fl9l) to be valid, we must have 
lim f(s) = +oo. In fact this is precisely what happens numerically: / 

satisfies equation (1201) to the left of s with initial condition /(0) = 0; / 
also satisfies (120]) to the right of s with condition lim f(s) = +oo. In 

s—>s 

particular, / does not depend on the size of the kink, only the fact that 
there is one. If there are multiple kinks, we get the same condition: / 
approaches infinity from the right side of the kink. 

Knowing a,i and 6j, we know M !J . Recall from ([HD the formula 

n _. m— 1 rn— 1 

m=maxij k=i 1=] 
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oo 



Pick any x and y in [0, 1] and s > max(x, y). For large n, choose % = nx, 

K 2 

j = ny_, and m = ns. Since a^ ~ 1 — ^\ except at the kink, we have 

n -1 " j 1 if s ^ (x, s) 

fc=nx [cos a if s e (x, s) 

for any x and s. On the other hand, we have 



ns-i /ns-i 

fc=nx \ fe=nx 



1 — _1 fc 

n ^-^ \n 
fc=nx 



O 




exp / f(a) da . 



Thus we have 
G(x,y) 

where 



(x, s)(f)(y, s)e~ f * f{a) da e- fi /(t) dT ds 



max(i,|/) 



(x,s) 



cos a x < s < s 

1 otherwise 



Thus we see G(x, y) > even if the curve has one or more kinks. So 
in the absence of gravity, the curvature is nonnegative in all sections. 

If there is gravity, then G(x, y) may be negative, if the string is 
pointing upwards at the fixed point. If we use a generalized curvature 
that incorporates the Hessian of the potential energy, then we will also 
get negative generalized curvature in this case. 
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